Nature Cancer

Follicular lymphoma (FL) is an incurable, prototypical relapse-remitting cancer, implying the existence of therapy-persistent cells that survive frontline treatment and seed disease recurrence1-3. However, these persister cells remain difficult to study directly in patients because immediate post-treatment sampling is ethically and practically challenging. Using a genetically defined mouse model that allows sampling of persistent cells immediately after frontline R-CHOP therapy, we prospectively isolate and functionally define relapse-founding cancer persister cells (CPC). The CPC is an IgM memory-like B-cell with high germinal center re-entry capacity. This state represents a discrete component of a heterogeneous residual pool indicating that residual disease is polytypic and that relapse potential may depend on which cells persist rather than on residual tumour burden alone. By integrating mouse CPC with human FL datasets, we show that an analogous transcriptional programme is detectable at diagnosis and is enriched in patients with inferior clinical outcome across independent cohorts4,5. These findings support the concept that relapse risk is linked to a conserved, genotype-agnostic CPC programme present before therapy. To explore therapeutic vulnerabilities, we developed a scalable in-vitro platform that models the CPC-like state and used it to identify sensitivity to histone deacetylase inhibition. Romidepsin and panobinostat killed CPC-like cells in-vitro, and decreased therapy-persistent cells after R-CHOP treatment in-vivo and in patient-derived organoids. Together, these data define a tractable CPC state in FL, with a validated clinical readout and an immediately testable therapeutic entry point, opening CPC-directed strategies for durable FL control.

Most high-dimensional studies of tumor-immune interactions focus on metastatic models, limiting insight into how immune remodeling in primary tumors shapes metastatic competence. Here, integrating single-cell RNA sequencing, CyTOF, and functional studies across metastatic (4T1) and non-invasive (EMT6) triple-negative breast cancer (TNBC) murine models, we define tumor state-specific immune programs that distinguish metastatic competence. Tumors with metastatic capacity uniquely drive early bone marrow expansion of CXCR2 neutrophils, which infiltrate primary tumors acquiring a CXCL2-producing phenotype that promotes EMT-associated cancer stem cell (CSC) plasticity. This program depends on TGF-{beta}/CEBPD-mediated induction of S100A9. Elevated CXCL2, together with G-CSF, establishes a feed-forward circuit that drives systemic neutrophil mobilization and recruitment to distant organs, where neutrophil-derived S100A8/A9 (calprotectin) promotes MET-driven CSC outgrowth and metastatic colonization. Clinically, gene signatures associated with CXCR2 neutrophils predict poor survival in TNBC patients, whereas monocyte/macrophage (CX3CR1) and T cell activation signatures correlate with improved outcomes. S100A9 ablation disrupts this cascade and enhances immunotherapy responsiveness, defining a TGF-{beta}/S100A9/CXCR2 axis linking immune remodeling, CSC plasticity and metastasis. HighlightsO_LIMetastatic TNBC engages a TGF-{beta}/C/EBP{delta}/S100A9 axis that expands CXCR2 neutrophils C_LIO_LINon-invasive EMT6 tumors retain a CX3CR1 monocyte/macrophage and T-cell landscape C_LIO_LICXCR2+ neutrophils in pre-metastatic niches suppress T cell response while promoting tumor cell proliferation C_LIO_LIS100A9 loss redirects myelopoiesis and potentiates anti-PD-L1 in TNBC models C_LI In BriefAlkan et al. dissect how tumor state programs the myeloid compartment in TNBC. Metastatic 4T1 tumors uniquely engage a TGF-{beta}/C/EBP{delta}/S100A9 axis driving CXCR2 neutrophil expansion and CXCL2/G-CSF-dependent systemic mobilization, coupling immune remodeling to EMT/MET cancer-stem-cell plasticity, while S100A9 loss restores CX3CR1 myeloid identity and unlocks checkpoint-inhibitor responsiveness.

Genomics-guided precision oncology has improved survival in cancer entities with actionable mutations but cannot capture oncogenic signaling that manifests at the protein level. Here, we report a prospective, real-world pan-cancer study profiling proteomes and phosphoproteomes of 1,998 tumor samples from adults and children with rare or advanced cancers enrolled in the German precision oncology programs DKFZ/NCT/DKTK MASTER, CATCH and INFORM and their molecular tumor boards (MTBs). We developed tumor proteome activity status (TOPAS) scores for 46 clinically relevant kinases, an immune activity score capturing antigen presentation and T-cell activation and identified therapeutically targetable cell-surface proteins for 94% of patients. These readouts enhance MTB recommendations by exposing actionable non-genomic kinase activity, refining interpretation of oncogenic genome alterations, and highlighting cell-surface treatment options. Three proof-of-concept analyses indicate clinical utility including kinase activity-stratified pazopanib response in sarcoma, immune activity score-tracked checkpoint-inhibitor outcomes pan-cancer, and a phosphoproteomic biomarker distinguishing EGFR-inhibitor response in chordoma.

Activating HER2 mutations are significantly enriched in both primary and metastatic invasive lobular breast cancer (ILC), with large public datasets of primary breast tumors linking them to a worse prognosis in ILC. Despite their oncogenic role, no FDA-approved therapies currently target HER2-mutant breast cancers. While the HER2-directed antibody-drug conjugate (ADC) trastuzumab deruxtecan (T-DXd) has shown efficacy in HER2-mutant non-small cell lung cancer, its activity in HER2-mutant ILC remains unknown. Using the Caris real-world database, one of the largest cohorts with survival data in advanced breast cancers, we confirmed that HER2 mutations are more prevalent in advanced ILC than in invasive breast cancer of no special type (NST) tumors, are associated with worse survival in both histologies, yet predict improved response to T-DXd across subtypes, highlighting the need for mutation-directed, histology-informed therapies. Using endogenous HER2-mutant ILC cell lines (UACC3133-S310F, BCK4-L755S) and CRISPR-engineered isogenic ILC models with clinically relevant HER2 mutations (S310F, V777L), we found these mutations drive HER2/HER3 hyperactivation and downstream signaling, conferring increased sensitivity to HER2 tyrosine kinase inhibitors (TKIs) and T-DXd. Mechanistically, HER2 mutants showed enhanced receptor ubiquitination, internalization, and lysosomal degradation upon T-DXd treatment, explaining the observed drug sensitivity. While combining T-DXd with neratinib or the HSP90 inhibitor ganetespib yielded synergistic effects in long-term growth assays, accompanied by increased HER2 ubiquitination, the concurrent hyperactivation of HER3 in HER2-mutant cells suggested that co-targeting HER3 could provide an effective alternative strategy. Accordingly, HER2-mutant ILC exhibited enhanced sensitivity to the HER3-directed ADC patritumab deruxtecan (P-DXd) or LJM716, a HER3-targeting antibody. We further uncovered a previously unrecognized mechanism of P-DXd beyond HER3 ligand blockade and payload delivery: P-DXd promotes HER2/HER3 association, increases HER2 ubiquitination, and enhances T-DXd internalization, resulting in potent synergy with T-DXd. Mechanistically, we identified HER3 extracellular domains I and II as essential for P-DXd binding and for mediating P-DXd-induced HER2/HER3 association, establishing a structural basis for this activity. In vivo, both T-DXd and P-DXd suppressed UACC3133 and BCK4 xenograft growth, with combination therapy trending toward greater efficacy and prevented regrowth of tumors. Extending these findings beyond HER2-mutant ILC, combination treatment with T-DXd and P-DXd demonstrated synergistic activity across multiple breast cancer models, including (i) HER2-amplified NST patient-derived organoids (PDOs) harboring hotspot HER2 mutations, (ii) HER2-wild-type NST PDOs with clinically intrinsic or acquired T-DXd resistance, and (iii) isogenic HER2-mutant ILC PDOs with experimentally induced resistance after prolonged T-DXd exposure. Collectively, these findings support HER2 as an actionable target in HER2-mutant ILC and position T-DXd-based regimens, particularly in combination with HER3 inhibition, as a promising therapeutic strategy for this underserved patient population.

Homologous recombination deficiency (HRD) is an actionable vulnerability found in a substantial fraction of human cancers, yet current HRD-directed therapies are limited by toxicity, incomplete responses, and acquired resistance. Many DNA-damaging agents were developed before DNA repair biomarkers were available, suggesting that abandoned agents may harbor previously unrecognized genotype-selective activity. Here, through a focused screen of DNA-damaging agents in isogenic homologous recombination-proficient and -deficient models, we identify CB1954, a decades-old nitrobenzamide aziridine prodrug, as highly selective for BRCA2-deficient tumor cells. CB1954 forms DNA interstrand crosslinks independent of HR status, but selectively induces DNA-damage signaling, apoptosis, and loss of clonogenic survival in HR-deficient cells. Targeted DDR CRISPR screening and isogenic validation define a distinct repair dependency for the Fanconi anemia and homologous recombination pathways, with limited dependence on mismatch repair or nucleotide excision repair. Genetic and pharmacologic perturbation of NQO2, the bioactivating enzyme for CB1954, reveals a bifurcated mechanism in which NQO2-dependent activation selectively contributes to HRD cytotoxicity, while aziridine-dependent lesions likely account for residual activity in HR-proficient cells. CB1954 exhibits favorable preclinical pharmacokinetic properties and genotype-dependent antitumor activity in BRCA2-deficient xenografts. These findings reposition CB1954 as a historically overlooked HRD-selective agent and demonstrate that biomarker-guided profiling of DNA-damaging agents can uncover new opportunities for precision oncology.

Rectal cancers are often treated with neoadjuvant chemoradiotherapy (CRT), yet 85% of patients do not achieve a pathological complete response. To identify the molecular determinants of CRT response, we profiled the single-cell signalling, DNA-damage, cell-cycle, apoptotic, and cell-fate responses of 2,769 patient-derived organoid cultures treated with CRT, cancer-associated fibroblasts (CAFs), and signal-rewiring agents. We find that CRT response is determined by stem cell-fate. CRT triggers comparable DNA-damage in isogenic proliferative (proCSC) and revival (revCSC) colonic stem cells, but proCSC retain damage and die whereas revCSC resolve damage and persist. Both CRT and CAFs drive proCSC to a common treatment-resistant revCSC fate and high revCSC predicts worse survival in patients. Pharmacologically constraining stem-cell plasticity increases CRT sensitivity, and Spatial Perturbation of ARrayed Tumour Assembloids (SPARTA) confirms YAP/TEAD inhibition improves chemotherapy responses in human stromal-tumour models. These results suggest that cancer cell-fate, not genotoxic damage itself, ultimately governs response to standard-of-care chemoradiotherapy. HIGHLIGHTSO_LIRectal cancer stem cell-fate determines chemoradiotherapy-induced apoptosis C_LIO_LIproCSCs retain DNA-damage and die, whereas revCSCs repair damage and persist C_LIO_LICAFs and chemoradiotherapy converge on a common chemo-radioresistant revCSC state C_LIO_LISPARTA reveals TEAD inhibition blocks DNA-repair persisters in stromal assembloids C_LI

Clinical-grade RAS inhibitors raise an unresolved question as to whether KRAS-alleles impose constraints on adaptive resistance that can be exploited therapeutically. Using daraxonrasib (RMC-6236), a multi-selective RAS(ON) inhibitor, we compared resistance mechanisms between KRASG12D and KRASG12R, alleles with fundamentally different RAS network dynamics. Daraxonrasib inhibited KRASMUT primarily through steric occlusion of effector binding, while engaging RASWT only modestly ([~]20%). KRASG12R is marked by its inability to transactivate RASWT, and it was observed that daraxonrasib resistant KRASG12R PDAC cells utilize EGFR/RASWT-GTP signaling as the dominant adaptive route. In contrast, KRASG12D resistance arose through retained KRASG12D-GTP signaling, with a decrease of cyclophilin A (CypA) protein, the binding partner required for daraxonrasib activity. The shift from KRASG12R dependence to the EGFR/RASWT conferred sensitivity to trametinib. We confirmed this clinically: a KRASG12R PDAC patient who progressed after 10 months on daraxonrasib showed intratumoral EGFR/RASWT activation, and rapid 3D-bioprinted patient-derived toroid modeling predicted sensitivity to trametinib-based combination therapy. Given the aggressive disease trajectory and lack of response to the two immediately preceding lines of therapy, sixth-line trametinib-based combination therapy achieved approximately 5 months of disease control. This patient ultimately achieved 40 months of overall survival, far exceeding the 8-12 month median for metastatic PDAC. Collectively, these data establish a framework in which allele-specific RAS network topology dictates the adaptive resistance landscape, enabling rational selection of targeted therapies with meaningful clinical benefit in metastatic PDAC. STATEMENT OF SIGNIFICANCEDaraxonrasib resistance mechanisms have allele-specific routes: CypA becomes downregulated in KRASG12D and reliance on EGFR/RASWT in KRASG12R. Rapid patient-derived toroids identified sixth-line targeted therapy strategies with an overall survival of 40 months.

Mucosal cannabidiol formulations are known regulators of the glioblastoma microenvironment, yet the underlying origin point triggering this stroma-remodeling efficacy remains entirely unknown. Here, by mapping innate cell trafficking pathways, we define a novel baseline neuro-immune-microbiome axis in orthotopic glioblastoma, characterized by diverse microbial communities, likely seeded via blood-brain barrier disruption, paired with dense infiltration of host mast cells and mature, crystalloid-containing eosinophils. Localized intranasal administration of a synthetic cannabidiol formulation achieved striking therapeutic efficacy, driving dramatic tumor regression. Mechanistically, high-throughput 16S rRNA sequencing and quantitative flow cytometry revealed this progression was subverted by taming the tumor ecosystem; cannabidiol restricted chaotic microbial diversity, selectively filtering the landscape toward Delftia and depleting Archaea, while simultaneously suppressing hyper-inflammatory host mast cell and eosinophil populations. This study builds upon established innate trafficking frameworks to present the first therapeutically targetable stromal-microbial axis in neuro-oncology.

Combination targeted therapy with BRAF/MEK inhibitors and immune therapy show promising therapeutic outcomes in melanoma; however, the development of drug resistance still represents a formidable challenge. Remaining unexplored is the possibility that BRAF/MEK inhibitors themselves inadvertently compromise the tumor immune microenvironment, limiting the efficacy of immunotherapy when it is used in combination with targeted inhibitors. Herein, we profiled the landscape of the BRAF regulatome identifying a novel transcription factor, TFAP2A, newly linking BRAF/MEK drug resistance to antitumor immunity. Specifically, we found that BRAF/MEK inhibitors significantly upregulate TFAP2A. Further, genetic disruption of TFAP2A overcomes BRAF/MEK-inhibitor resistance, promotes stromal enrichment, and enhances intratumoral infiltration of macrophages in an immune-compromised mouse model. In a syngeneic mouse model, TFAP2a knockout not only suppresses tumor growth but also induces potent anti-tumor tertiary lymphoid structures (TLSs). Single cell transcriptomics revealed that the absence of TFAP2A shapes the antitumor microenvironment with an influx of M1-like macrophages, CD8+ T cells and mature dendritic cells. By identifying TFAP2A as a shared driver of both targeted therapy resistance and immunosuppression, our work offers a one-stone-two-bird strategy to overcome drug resistance and elicit antitumor immunity.

Macrophages are innate immune cells of embryonic or adult origin with tissue specific roles in homeostasis, disease surveillance, and wound repair that can be co-opted to promote tumor growth and spread1-11. An understanding of the specific roles of macrophage subsets in lung tumor initiation and progression could promote new therapeutic approaches for this deadly disease. Here, we show that KRASG12D mutations in lung epithelium drive proliferation of resident, embryonically-derived alveolar macrophages, which then promote tumor cell proliferation and protection from ferroptosis, leading to tumor progression. Using genetically engineered mouse models of mutant KRASG12D non-small cell lung cancer12,13, we found that alveolar macrophages accumulate by proliferation in response to tumor cell-secreted IL-34, recapitulating events observed in late embryonic lung development. Tumor alveolar macrophages in turn drive IGF-1-dependent tumor cell proliferation. Neutralization or deletion of IL-34 suppresses IGF-1 expression, reduces macrophage and tumor cell proliferation and inhibits tumor progression. High IL34 and IGF1 correlate with poor survival in KRASG12D/V lung adenocarcinomas and in other solid tumors, indicating that bi-directional proliferative signaling between resident macrophages and tumor cells can drive human lung tumor progression. These studies identify resident macrophage-tumor cell interactions as key interception points for lung cancer therapy.

Hematologic malignancies are diagnosed through a fragmented, sequential workup of morphology, immunophenotyping, cytogenetics, and molecular testing that can take days to weeks and is unavailable at many centers. DNA methylation profiling has transformed central nervous system tumor diagnosis, yet hematologic classifiers have remained confined to narrow acute leukemia panels. Here we present Lamprey, a deep-learning methylation classifier spanning 86 hematologic malignancy entities, trained on a reference cohort of 8,544 patients and deployed directly from nanopore sequencing. A depth-aware training framework allows confident classification from the first minutes of a run. Against blinded integrated reference diagnoses across retrospective, external, and prospective cohorts, Lamprey exceeded 98% accuracy among classified cases. Lamprey reaches a confident call within minutes, and cost as little as $82 per sample. Lamprey consolidates a sequential diagnostic workup into a single, rapid, same-day molecular readout.

Hypermutational processes, including those driven by the APOBEC3 family of cytidine deaminases, generate abundant neoantigens yet give rise to tumors that evade immune recognition. Here, using multi-omics analyses followed by functional validation, we identified a tumor-intrinsic immune-cloaking mechanism in neoantigen-rich epithelial cancers, characterized by coordinated suppression of antigen presentation, immune-recruiting cytokines and immune-checkpoint programs. In bladder cancer, genome-wide copy-number analysis identified recurrent 1q23.3 amplification as a genomic feature of a neoantigen-high/CD8-low tumor state. Within this locus, NECTIN4 emerged as the dominant candidate effector, outperforming extrachromosomal DNA status as a predictor of immune-neoantigen discordance. Similar associations were observed across breast and lung cancers. Functional studies demonstrated that NECTIN4 was sufficient to establish a T-cell-poor tumor microenvironment and confer resistance to PD-1 blockade in immunocompetent mice. Mechanistically, NECTIN4 engaged a DDR1-SHP2 axis that suppressed STAT1 phosphorylation, silencing tumor-cell immune-engagement programs. NECTIN4 blockade restored STAT1 activity and reduced tumor growth, indicating that the cloaked state is pharmacologically reversible. Mutational signature, breakpoint motif, timing and clonality analyses, together with APOBEC3B expression and germline genetic evidence, linked APOBEC3-mediated mutagenesis to recurrent 1q23.3 amplification encompassing NECTIN4. These findings reveal how neoantigen-generating mutational processes can be coupled to structural genome evolution to enable tumor-intrinsic immune cloaking through a therapeutically targetable NECTIN4-DDR1-SHP2 axis.

GD2-CAR T cell therapy has demonstrated clinical benefit in patients with H3K27M+ diffuse midline glioma (DMG), but the durability of response has been limited in many patients1,2. To identify mechanisms of therapeutic resistance, we conducted longitudinal single-cell RNA and TCR sequencing of cerebrospinal fluid (CSF) lymphocytes from DMG patients receiving intravenous followed by sequential intracerebral GD2-CAR therapy, with lymphodepleting chemotherapy administered once prior to the start of CAR T cell therapy (NCT04196413). CSF GD2-CAR T cells manifested limited persistence and clonal expansion, while non-engineered CSF lymphocytes underwent significant clonal expansion and repertoire stabilization, ultimately dominating the CSF immune compartment. Concurrently, peripheral blood CD4+ and CD8+ T cells manifested anti-CAR immune reactivity targeting epitopes enriched within murine-derived and engineered junctional regions of the CAR construct. This was associated with appearance of circulating Human Anti-CAR Antibodies (HACAs) that bound cells expressing the GD2-CAR, as well as clonal expansion of CSF B cells which produced HACA which impeded the cytotoxic activity of GD2-CAR T cells. In several cases, appearance of circulating HACA temporally correlated with disease progression and across the patient population, and levels of circulating HACA inversely correlated with circulating CAR T cell persistence. These findings reveal robust induction of systemic and CNS adaptive T cell and B cell responses to GD2-CAR T cells following intravenous then sequential intracerebroventricular GD2-CAR therapy and provide strong evidence that anti-CAR immunity is a significant contributor to therapeutic resistance in this setting.

Triple-negative breast cancer (TNBC) exhibits pronounced molecular heterogeneity, yet the majority of published transcriptomic prognostic signatures suffer from limited reproducibility and have not achieved clinical translation. We systematically benchmarked 62 published TNBC prognostic signatures across 6 independent cohorts (n=1,357) using a unified analytical framework spanning multiple scoring algorithms, survival endpoints, and threshold strategies. While 17 signatures demonstrated consistent univariate prognostic associations, only 4 remained independently prognostic after adjustment for clinicopathological variables, and none achieved robustness across all analytical conditions-underscoring the fragility of existing classifiers. Leveraging genes with concordant survival associations across all 6 discovery cohorts, we identified a reproducible 15-gene directionally concordant gene set (DCGS) signature and distilled it into MetaSig-EFS, a 13-gene prognostic model optimized using a cohort-aware DeepSurv framework. MetaSig-EFS demonstrated robust cross-cohort generalizability, achieving validation concordance indices of 0.89 in GSE19615 and 0.69 in JBordet, with corresponding 3- and 5-year time-dependent AUROCs of 0.89 and 0.96 in GSE19615 and 0.80 and 0.71 in JBordet, respectively, while retaining independent prognostic value across established TNBC molecular sub-typing systems. Leveraging the TAHOE-100M transcriptomic perturbation atlas, we systematically prioritized candidate therapeutics through large-scale drug repurposing, identifying Paclitaxel as the top-ranked compound, followed by Venetoclax and Tucatinib. Together, these findings provide biologically informed and clinically actionable therapeutic hypotheses that extend the translational utility of our validated prognostic framework for TNBC risk stratification.

BackgroundGlioblastoma (GBM) invasion is clinically decisive but difficult to model systematically. Existing patient-derived xenograft (PDX) resources rarely couple reproducible in vivo invasion phenotypes with matched multi-omic profiles at scale, limiting mechanistic insight and phenotype-informed therapeutic hypotheses. MethodsWe established the HGCC Phenobank, comprising 65 patient-derived GBM stem-like cultures with matched multi-omic profiling and orthotopic engraftment in 449 mice. Blinded histopathology quantified ten invasion traits per case. These phenotypes were integrated with RNA sequencing, DNA methylation, and mass-spectrometry-based proteomics. Multi-Omic Factor Analysis (MOFA) identified latent molecular programs. Phenotype-specific RNA signatures were matched to LINCS drug-perturbation profiles and validated in 3D gliomasphere and ex vivo brain-slice assays. ResultsTwo dominant, reproducible invasion modes emerged across models: diffuse parenchymal infiltration and perivascular/condensed growth. Proneural cultures formed more aggressive tumors in immunodeficient mice, and mouse survival showed a modest correlation with patient survival in matched cases (Pearson r = 0.1832, p = 0.045). MOFA identified 15 latent factors; Factor 1, enriched for ASCL1/OLIG1/OLIG2 programs and associated with TP53/DCHS2/WNK2 alterations, was linked to increased tumor formation, diffuse invasion, and shorter mouse survival, and stratified GBM patients in TCGA and in our matched patient cohort. Drug-signature matching separated mechanisms targeting diffuse versus perivascular invasion. Experimental validation confirmed phenotype-selective sensitivities, and inhibitors PIK-75 and buparlisib suppressed invasion dynamics across representative models in 3D and brain-slice assays. ConclusionsThe HGCC Phenobank provides the first openly available PDX resource that systematically links GBM invasion phenotypes to multi-omic programs and therapeutic predictions. This framework enables reproducible model selection, mechanistic dissection of invasion modes, and phenotype-guided therapeutic discovery. Key PointsO_LIDiffuse and perivascular invasion define orthogonal GBM axes C_LIO_LIASCL1/OLIG factor links initiation, diffuse growth, and survival C_LIO_LIPhenotype-matched drugs validated; PIK-75 and buparlisib curb invasion dynamics C_LI Importance of the StudyGlioblastoma invasion varies substantially between patients, yet existing patient-derived xeno-graft resources rarely combine reproducible in vivo phenotyping with matched multi-omic profiling at scale. The HGCC Phenobank addresses this gap with standardized, blinded scoring of ten invasion traits across 449 orthotopic xenografts from 65 molecularly characterized GBM stem-like cultures, integrated with transcriptomic, methylomic, and proteomic data. We identify two dominant, reproducible invasion modes and a cross-modal neurodevelopmental program, the ASCL1/OLIG1/2-associated Factor 1, that links tumor initiation, diffuse growth, and survival in mice, and stratifies GBM patients in TCGA and in our matched patient cohort. In a spatially resolved xenograft section, Factor 1 signal localizes to the invasive tumor periphery. By matching phenotype-specific RNA signatures to drug-induced transcriptional responses, we show that invasion phenotypes nominate selective vulnerabilities, exemplified by PIK-75. This openly shared resource enables reproducible model selection, mechanistic dissection of invasion programs, and phenotype-guided therapeutic discovery.

Cancers with similar histology often exhibit divergent clinical behavior, reflecting molecular heterogeneity not captured by current classification systems. Although driver mutations are central to tumorigenesis, their broader systems-level consequences have not been systematically leveraged. We integrated genomic and transcriptomic data across cancers to define driver mutation signatures (DMS), coordinated transcriptional programs associated with cancer driver mutations. From 121 candidate drivers, we derived 90 robust signatures and quantified their activity in individual tumors using mutation signature scores (MSS). DMS analysis revealed a hierarchical organization of tumors into molecular subgroups that transcended tissue boundaries while preserving driver-associated features. Continuous MSS profiles further defined high-resolution molecular fingerprints for individual tumors. DMS provides a quantitative framework for tumor classification and patient stratification and links driver-associated programs to potential therapeutic vulnerabilities. Together, these findings establish a pan-cancer molecular taxonomy that bridges genotype and phenotype and may inform precision oncology.

Somatic copy-number amplification is pervasive in cancer, and the genes it carries are candidate drug targets - but only those whose amplification is transmitted to accessible surface protein can be reached by an antibody-drug conjugate (ADC). We build an integrated map of copy-number-to-protein transmission across six tumour types and ask, for every amplified gene, whether its dosage reaches the surface. Copy number transmits to mRNA (median per-gene r = 0.21) but is attenuated at the protein level in 85% of genes, and the mRNA ranking is largely preserved to protein (rho = 0.70); the ranking is set principally at the chromatin/transcription step - among directly measured regulatory inputs, promoter DNA methylation and tumour chromatin accessibility each explain about an order of magnitude more of the transmission variance than gene structure, and do so complementarily. Critically, transmissibility is a stable, gene-intrinsic property: it is predictable from gene properties alone, with no proteomic input, at a leave-gene-out rank correlation of 0.52 (R2 = 0.29); it is not positional (holding out whole chromosome arms changes accuracy by 0.001); and it transfers across lineages (Kendall W = 0.97 across leave-one-lineage-out refits). This licenses a predictor that nominates surface targets in cancer types that lack a tissue-referenced proteome, combining direct protein measurement where it is available with prediction where it is not. Requiring co-elevation on a recurrent amplicon with measured transmissibility and an accessible extracellular ectodomain nominates 22 surface antigens on 18 distinct recurrent amplicons across four cancer types (renal, endometrial and both lung subtypes) - for example ITGB8+TSPAN13+TTYH3 on lung 7p, NCSTN+HSD17B7+MPZL1 on 1q (recurrent in several types), the transferrin receptor TFRC on squamous 3q, and FZD1 on clear-cell renal 7q; 21 of the 22 are non-driver passengers and 10 are confirmed on the experimental Cell Surface Protein Atlas. In single malignant cells, against a null that controls for per-cell sequencing depth, the co-detected constructs sit at a modest 1.05-1.45x above independence (p < 0.001, donor-block bootstrap intervals clear of 1.0), and at binding-relevant thresholds the normal-tissue co-expression collapses - so an avidity AND-gate that binds stably only where the antigens co-occur would spare normal cells that carry only one. Observed transmissibility itself transfers strongly between the two lung subtypes ({rho} = 0.88) and remains positive across distant lineages, consistent with the shared cell-of-origin regulation the map implies. Single-cell co-detection is demonstrated wherever a malignant single-cell atlas exists (both lung subtypes and glioblastoma - the latter entirely from prediction, using no GBM surface-abundance measurement); the remaining cohorts are nominated on the same genetic and topological evidence. The result is a pan-cancer, confidence-tiered catalogue of multi-antigen ADC co-target sets with a concrete plan to test them.

Targeting oxidative phosphorylation (OXPHOS) represents an attractive therapeutic strategy in acute myeloid leukemia, which exhibits exceptional dependence on mitochondrial respiration compared to normal hematopoietic cells. However, clinical attempts to exploit this vulnerability have been limited by on-target toxicity to healthy tissue. Here, we comprehensively compare the cellular consequences of inhibiting distinct nodes of the electron transport chain in AML. We demonstrate that selective inhibition of the F1 subunit of ATP synthase with EB2023 (ammocidin A) delivers an energetic stress to AML cells without the profound redox stress that characterizes complex I inhibition, preventing NAD/NADH imbalance and allowing continued TCA cycling. Further, the duration of OXPHOS inhibition is transient in nature in vivo, a finding revealed through pharmacokinetic and serial pharmacodynamic monitoring of AMPK phosphorylation accompanied by OPA1-mediated mitochondrial structural remodeling that primes AML cells for BCL2 inhibitor synergy. EB2023 in combination with venetoclax demonstrates potent anti-AML activity across cell lines and patient-derived xenograft models at doses that spare normal hematopoietic progenitors and avoid the neuropathy and sustained detrimental systemic metabolic rewiring in healthy tissues associated with prior efforts to target OXPHOS. These findings establish F1-selective ATP synthase inhibition as a clinically actionable therapeutic strategy in AML and establish the duration of OXPHOS inhibition as a critical and previously underappreciated determinant of therapeutic index.

Clear cell renal cell carcinoma (RCC) patients have multiple approved therapies, including anti-angiogenesis tyrosine kinase inhibitors (TKIs) and immuno-oncology therapies (IO), but lack a clinically validated biomarker. Using RNAseq data from the IMmotion 151 clinical trial (IM151)1-3, 7 RCC biologic clusters have been defined3. Several groups have attempted to predict these clusters on various RCC datasets4,5 with negative results, suggesting that the association with therapy response from IM151 could not be reproduced. We hypothesized that the specific approaches used to generate cluster predictions led to the misinterpretation of findings. Both published models used standardization (z-scores) to normalize the data within their patient cohorts, imposing an expected gene expression distribution in which [~]50% of patients have apparently higher-than-average expression, artificially impacting the proportion of cluster assignments and leading to potential misclassification. We developed a machine learning (ML) model, IRIS-RCC (Individual RNA-seq Intrinsic Subtyping for RCC), to predict treatment (TKI vs. IO) for patients using an individual-sample model using the IM151 trial (N=823) and validated the model on the JAVELIN Renal 101 trial (JR101; N=726)3,6. Our method normalizes gene expression within a given sample using ratiometric expression. This method results in different cluster assignments for individual tumors and distinct clinical correlations. An additional advantage of individual-sample predictions is that they can be readily applied in a prospective setting, where patients must be classified one at a time. IRIS-RCC is currently being validated in a prospective biomarker-driven Phase II clinical trial (OPTIC RCC).

Hormone receptor-positive, HER2-negative breast cancers are often poorly infiltrated by immune cells and derive limited benefit from current immunotherapy strategies. Here, using paired tumour samples from the randomised window-of-opportunity D-BIOMARK trial (NCT03691311), we investigated the immunomodulatory effects of denosumab in early luminal breast cancer. Short-term preoperative denosumab reduced tumour-cell proliferative transcriptional programs and immunosuppressive features of the local tumour microenvironment, enhancing innate and adaptive immune activation and altering circulating cytokine profiles. High-resolution spatial transcriptomics revealed coordinated remodelling of tumour, immune, fibroblast and endothelial compartments after treatment. Denosumab reduced immune-tumour spatial separation and enhanced T cell activation, accompanied by a shift from matrix-associated tumour programs towards increased tumour-T cell communication. Copy number-informed tumour-state inference further identified a reduced representation of genomically complex, immune-poor tumour subclones after treatment. Together, these findings identify denosumab as a modulator of tumour-microenvironment crosstalk and support RANKL blockade as a strategy to render immune-poor luminal breast tumours more permissive to immune engagement.